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ABSTRACT. Generalized Standard Materials are governed by maximal cycli- 
cally monotone operators and modeled by convex potentials. Gery de Saxce's 
Implicit Standard Materials are modeled by biconvex bipotentials. We ana- 
lyze the intermediate class of n-monotone materials governed by maximal n- 
monotone operators and modeled by Fitzpatrick's functions. Revisiting the model 
of elastic materials initiated by Robert Hooke, and insisting on the linearity, 
coaxiality and monotonicity properties of the constitutive law, we illustrate that 
Fitzpatrick's representation of n-monotone operators coming from convex anal- 
ysis provides a constructive method to discover the best bipotential modeling a 
n-monotone material. Giving up the symmetry of the linear constitutive laws, 
we find out that n-monotonicity is a relevant criterion for the materials charac- 
terization and classification. 



1. Introduction 

Standard Materials are modeled by differentiable potentials. Mainly to capture 
set-valued constitutive laws, for example plastic flow rules ( l24l 1401 ) they were 
extended to the so-called "Generalized Standard Materials" (GSM) modeled by 
lower semi-continuous (lsc) convex potentials ([28, 29 1). However, this extension 
failed to describe Coulomb's dry friction law. 

In 1991, considering such an implicit constitutive law, Gery de Saxce and Zhi- 
Qiang Feng ( 11471 148 49 1) proposed a new generalization which they called "Im- 
plicit Standard Material" (ISM). This new class of materials is modeled by a point- 
to-point function which they called a bipotential. In the particular case of a GSM, 
the bipotential becomes simply the sum of the potential and its conjugate. 

For a given GSM, a theorem due to R. Tyrrel Rockafellar ([44 46 1) and Jean- 
Jacques Moreau ( BP ) provides a constructive method to retrieve the potential 
from its subdifferential. A similar, crucial question for a given ISM is: how to 
retrieve the bipotential from the implicit constitutive law? 

Independently, in 1988, Simon Fitzpatrick ([25 1), in order to simplify the study 
of monotone operators (HE)), made a proposal to replace these multifunctions 
by point-to-point functions, nowadays called Fitzpatrick's functions. 
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Is there a relation between Gery de Saxce's bipotentials representing ISM con- 
stitutive laws and Fitzpatrick's functions representing maximal monotone multi- 
functions? Can this last representation coming from convex analysis provide a 
method to construct the best bipotential modeling a given ISM? The aim of our 
paper is to discover the largest class of ISM for which the answers to both ques- 
tions are positive. 

Revisiting the modeling of elastic materials initiated by Robert Hooke, we il- 
lustrate Fitzpatrick's method on the example of linear coaxial constitutive laws. 

2. Standard Materials 

A constitutive law relating a strain-like variable x belonging to a real Banach 
space X (with norm || • ||) and a stress-like variable y belonging to the topological 
dual space Y = X* (with the duality product (•, •)), is a subset of the product set 
X x Y. This subset can be regarded as the graph 

G(T) = {(x,y)eXxY \yeTx} 

of a multivalued operator T : X — > 2 Y . In finite dimension, when this subset is 
a maximal lagrangian submanifold of the linear space X x Y (made symplectic 
by the canonical Darboux 2-form), there exists a differentiable function cf>, called 
"potential", such that the constitutive law reads 

y = grad</>(x). 

If, additionally, this potential is convex, the inverse constitutive law reads 

x = gr&d (f)*(y) 

with cj)* the Legendre transform of the potential 4> (called "conjugate potential"). A 
material whose behavior can be described by a differentiable potential is referred 
to as a "Standard Material" (SM). 

3. Generalized Standard Materials 

For many materials, the relation between x and y is multivalued. Giving up 
the differentiability of the potential <j>, but retaining its convexity and its lower 
semi-continuity, a large class of materials, called "Generalized Standard Materi- 
als" (GSM), can be described by one of the following three equivalent constitutive 
laws: 

(i) y G d<j>(x) 

(ii) x g d<f>*(y) 

(Hi) <j)(x) + <j)*(y) = (x,y). 

Remark 3.1. The convexity of the potential cf> allows to express the conjugate 
potential <j>* as a supremum ( ATI ) 

<p*(y) = sup [(x,y) - (j)(x)] . 
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Remark 3.2. The subdifferentials 



d<j>(x) 



{y G Y | V£ G X, > <f>(x) + {£-x, y)} 



and 



d<fr*(y) 



{xeX\V V e Y, & (77) > <t>* (y) + (x,r]- y)} 



generalize (ET)) the gradients of the potentials <f> and cf>* when these convex po- 
tentials are not differentiable. Their elements are called subgradients. At a point 
x (respectively y) where <j> (respectively <fi*) is both convex and differentiable, the 
set of subgradients d<j)(x) (respectively d(f>* (y)) reduces to the unique gradient at 
x (respectively at y). 



4.1. Bipotentials. The equality {Hi) of Section[3]can be regarded as an extremal 
case of Fenchel's inequality 



For modeling the dry friction phenomenon or the behavior of materials such as 
clays, Gery de Saxce noticed that it was fruitful to weaken this inequality to 



Foregoing to separate the function b(x, y) into a sum <fi(x) + <f>* (y) of two poten- 
tials, he called it a bipotential. The bipotentials b(x,y) are assumed to be 

(i) convex and lsc in x 

(ii) convex and lsc in y 

(iii) bounded from below by the duality product: b(x, y) > (x, y) 
as it is the case for separable bipotentials. 

4.2. Implicit Standard Materials. A material whose behavior can be described 
( lfT9ll20ll2Tl ) by one of the following three implicit constitutive laws 

(iv) y belongs to the subdifferential of &(£, y) with respect to £ at x 

(v) x belongs to the subdifferential of b(x, -q) with respect to r\ at y 



4. Implicit Standard Materials 



(f>(x) + <f>*(y) > (x,y). 



(vi) b{x,y) = (x,y) 
is referred to as an "Implicit Standard Material" (ISM). 
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4.3. Examples of Implicit Standard Materials. The Implicit Standard Mate- 
rial model proved to be relevant for describing many non-associated phenomena 
([50] and references contained therein): unilateral contact with Coulomb type dry 
friction ( I32ll33ll47ll49l ). generalized Dmcker-Prager plasticity (|50|), modified 
Cam-Clay model (031 1501 15111 ). non-associated plasticity of soils (10 ED), 
non-linear kinematical hardening rule for cyclic plasticity of metals ([1 36 1), 
Lemaitre's plastic-ductile damage law ( l35l ). and shakedown analysis of non- 
standard elastoplastic materials ( irT6l[T7l[T8ll5l"l ). 

In the next sections ([5] 17] [8]) we recall some notions and results on monotone 
multifunctions (also called multivalued operators or set-valued maps) that are well 
suited to model the behavior of non-associated materials. 

5. Basic facts on monotone constitutive laws 

5.1. Monotonicity. 

Definition 5.1. A constitutive law associated to a multifunction T is monotone 

(im[I2IH]|45]|59))if 

[yi G Txi and y 2 G Tx 2 ] [(x 2 ~ x 1 ,y 2 ~ y\) > 0]. 
It is strictly monotone ( ||59l ) if, additionally, 

[yi G Txi, y 2 G Tx 2 and x 2 ^ x%] [(^ 2 - xi,V2 - Vx) > 0]. 
Example 5.2. If A is a positive linear mapping from X to Y, i.e. 

Vx G X, (x, Ax) > 0, 
then the single-valued multifunction T defined by Tx = {Ax} is monotone. 

Example 5.3. If T is the subdifferential of a convex lsc potential (f>, then T = dcf> 
is monotone (ll8l l4Tl0 . However, a monotone multifunction is not necessarily the 
subdifferential of a lsc convex potential. 

5.2. Maximality. 

Definition 5.4. A monotone multifunction T is maximal if there is no monotone 
proper enlargement of T (||8l[TT1[T2ll4Tll59l). 

Remark 5.5. Maximal monotone multifunctions are monotone multifunctions 
with graphs that cannot be enlarged without destroying monotonicity. To estab- 
lish that a monotone multifunction is maximal, one must prove that 

(x,y) # G{T) =► [3( Xl , yi ) G G(T), (x -x x ,y- Vl ) < 0]. 

Remark 5.6. The maximality assumption is equivalent ([42, 52]) to the statements 

[(x, y) G X x Y and V(xi , y x ) G G(T), (x — x\,y — y\) > 0] => y G Tx, 
[(x, y) G X x Y and inf (x - x u y - yi ) > 0] => {x, y) G G(T). 

Lemma 5.7. If T is a maximal monotone multifunction, then 
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(i) (x, y) e G(T) =*> M yieTxi (x - x x , y - y x ) = 
(ii) (x, y) g G(T) =*> M yieTxi (x- xi,y-yi) < 

(Hi) yx e X, Vy e Y, ini VieTxi {x - x x ,y - yi) < 0. 

Proof. This lemma is a quantitative version ([42, 52 1) of Remark l5.6l Assertion 
(i) follows by taking (x x ,yx) — (x, y), assertion (ii) is immediate from the defi- 
nition of maximal monotonicity, and assertion (Hi) follows from (i) and (ii). □ 

The following example illustrates the important class of linear maximal mono- 
tone multifunctions. 

Example 5.8. Let X and Y be two instances of the same Hilbert space. If A 
is a positive linear (generally not symmetric) mapping from X to Y (Example 
15.2b . then the (single valued) monotone multifunction T defined by Tx — {Ax} 
is automatically maximal monotone (|6 42 52|). 

Proof. To show this, we must prove that 

z = y - Ax ^ =>■ [3xi £ X, {x - x x ,y - Ax x ) < 0]. 

Define x x — x + Xz where A is a scalar to be chosen later; then, by linearity, 

y — Ax\ = z — XAz, and 

(x — x\, y — Axi) — A(A(z, Az) — (z, z)). 
If (z, Az) = 0, the expression above is negative for any A positive. If (z, Az) is 

(z z) 

positive, the expression above is negative for < A < A'^v ■ In this second case, 
we can notice that according to Lemma l5l7l (ii) 

1 (z z) 2 

inf (x-x 1 ,y-y 1 ) < inf(A 2 (z, Az) - \{z, z)) = ~~;Ta \ < °- 

yiGT^i A 4 [Z, AZ) 

We recognize the inverse of Rayleigh's quotient 'f'v of the symmetric part S of 
A. Actually, if S is positive definite, then (Proposition 19.2b 

inf (x-xi,y-yi) =--(z,S~ 1 z) <Q. 

yiGTx! 4 

□ 

Thus, in the special case of linear mappings, the monotonicity implies the maxi- 
mal monotonicity. Surprisingly, the continuity is also automatically ensured ([52 |). 
However, it is worth mentioning that the situation is not so simple when A is only 
defined on a proper linear subspace of X r il42ll52l ). 
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5.3. Fitzpatrick's function. 

Definition 5.9. Let T be a maximal monotone multifunction. The associated Fitz- 
patrick function Ft,2 is defined (J25)) by 

Fr,2(x,y) = (x, y) - inf (x - x±, y - yi). 

Proposition 5.10. Fitzpatrick's function is bounded from below by the duality 
product 

F T ,i{x,y) > (x,y) 
and equality is attained if and only if y G Tx. 

Proof. It is a direct application of Lemma \5J\ □ 

Theorem 5.11. The Fitzpat rick function Ft, 2 represents the maximal monotone 
multifunction T: 

(i) G(T) = {(x,y) eXxY\ F T ,2{x,y) = (x,y)} 

(ii) outside G(T), FT y 2{x,y) > (x,y) . 

Proof. It is a recast of Proposition l5.10l □ 

Remark 5.12. Ft ,2 is globally lsc and convex. 
Proof. From the definition 

F T ,2{x,y)= sup [(x,yi) + (xi,y) - (xi,yi)] 

(x uV1 )eG(T) 

is the supremum of a family of continuous affine real-valued functions and there- 
fore is convex and lsc on the product space X x Y. □ 

Remark 5.13. With the duality product 

(((xi,yi), (v2,x 2 ))) = (xi,y 2 ) + (x2,yi) 

between X xY and Y x X, Fitzpatrick's function can be regarded as a Legendre- 
Moreau-Rockafellar transform. 

Proof. From the definition 

F T ,2{x,y)= sup [(((£1,2/1), (y, x))) - (xi,Vi)] 
(n,vi)eG{r) 

= sup [{{{xi,yi),{y,x)))-{x 1 ,y 1 )-i G (T)(x 1 ,y 1 )] 

which is nothing else than the value at (y, x) of the conjugate of the function 
(., .) + Ig(t)> where i'g(t) is the indicator function of the graph G(T). □ 

Example 5.14. Let X and Y be two instances of the same Hilbert space. Let 
S be a symmetric linear mapping from X to Y. If S is positive semi-definite, 
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then Fitzpatrick's function of the (single valued) multifunction T defined by Tx 

{Sx} is 



Ft,2( x >v) 



±(y-Sx,$ if (y - Sx) e R(S) 
and S£ = \{y - Sx) 
if (y - Sx) £ R(S) 



where R(S) is the range of S. 



Proof. The supremum with respect to x\ of the concave function (x, y) + (x\ — 
x, y — Sxi) is attained when y — Sx% = S(x% — x) or y — Sx = 2S(x\ — x). 
The result follows with £ = x\ — x as any solution of the linear equation 5£ = 
\ (y — Sx) when y — Sx belongs to the range R(S) of S. □ 

6. Basic facts on ?t,-monotone constitutive laws 
6.1. n-monotonicity. 

Definition 6.1. For n > 2, a constitutive law associated to a multifunction T is 
n-monotone if ([4]) 



22(x i+ i -Xi,yi) < 



for any n elements {x t ,yi) of G(T) and (x n+1 ,y n+1 ) = {xi,yi). 

Remark 6.2. We note that 2-monotonicity reduces to ordinary monotonicity. 

Remark 6.3. For n > 2, (n + l)-monotonicity implies n-monotonicity. 

Proof. Restrain the choice of the (n + 1) pairs in the graph G(T) by choosing 
(x n+ i,y n+ i) = (x n ,y n ) and (x„ +2 , y„ +2 ) = (xi,yi). □ 



where e is a nonzero 



Example 6.4. A 2-monotone example and a 3-monotone counterexample. 

\ 1 — 2el 

Let X = Y = R 2 and A be the 2 x 2 matrix _ 

2e 

scalar, then the (single valued) multifunction T defined by Tx = {Ax} is 2-mono- 
tone but not 3-monotone. 



Proof. The symmetric part S 



1 





of A is positive semi-definite, therefore T 



is 2-monotone but not strictly 2-monotone. Take x\ 



s 




"0" 










, X2 — 


1 


, x 3 






then yi = Aa;i = e 



2/2 = Ax 2 = -2e 



"1" 










. 2/3 






(x 3 - x 2 ,y 2 ) + (xi - x 3 ,y 3 ) reduces to 



— s 




"l" 


1 




2e 



, and (x 2 - xi,yi) + 
e 2 which is positive. 



With a positive semi-definite symmetric part and a nonzero skew symmetric 
part, the linear mapping A cannot generate a 3-monotone multifunction T. □ 
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6.2. n-monotonicity of positive semi-definite symmetric linear mappings. 

Proposition 6.5. Let X and Y be two instances of the same Hilbert space, and 
S a symmetric linear mapping from X to Y. If S is positive semi-definite, then 
the (single valued) multifunction T defined by Tx — {Sx} is n-monotone for any 
n > 2. 

Proof. For n = 3, in the inequality of the Definition 16.11 let us regard x 3 as the 
origin in X and set Xi — x 3 + Zi. This translation leads to force 

(zi, Szi) + (z2, Sz 2 ) - (z 2 , Szi) 

to be non-negative. Necessarily, S has to be positive semi-definite (take z 2 = 0) 
and we can extract its positive semi-definite square root S*. There is no other 
condition for the 3-monotonicity because 

(shx, si Zl ) + (sh 2 , sh 2 ) - (sh 2 , sh^ 

is non-negative as well as the 2x2 matrix 

" 2 -1" 
-1 2 ' 

For n = 4, in the inequality of the Definition l6.il let us regard X4 as the origin in 
X and set xi = X4 + Z{. This translation leads to force the non-negativity of 

(zi, Szi) + (z 2l Sz 2 ) + (z 3 , Sz 3 ) - (z 2 , Szi) - (z 3 , Sz 2 ). 

Necessarily, S has to be positive semi-definite (take z 2 = and z 3 = 0). There 
is no other condition for the 4-monotonicity because the above expression is non- 
negative as well as the 3x3 matrix 

"2 -1 " 
-1 2 -1 . 
0-12 

With the same arguments, one can prove that T is n-monotone for every n > 2. 
This result is consistent with the fact that T is the subdifferential of the convex 
potential <fi defined by <fi(x) = h{x, Sx) (see Example [jO} . □ 

6.3. 3-monotonicity of non symmetric linear mappings. 

Proposition 6.6. Let X = Y be a Hilbert space, and A a linear mapping from 
X to Y. A is not symmetric, the skew-symmetric part W of A is nonzero, and 
the symmetric part S of A is assumed to be positive definite. The ( single valued) 
multifunction T defined by Tx — { Ax} is strictly 2— monotone (Example 15.21 ). 
Furthermore, if the linear mapping 

H 3 = S- -A*S~ 1 A = -S - ^WS^W* 
4 4 4 

is positive, then the multifunction T reveals to be 3-monotone. 
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Proof. For n = 3, in the inequality of Definition l6.il let us regard x$ as the origin 
in X and set Xi = x$ + z;. This translation leads to force the non-negativity of 

2(zi, Szi) + 2(z 2 , Sz 2 ) - (z 2 ,Azi) - (A*z 2 ,z 1 ). 

Let us set z 2 — Az±, with A a linear mapping from X to X. Then, the non- 
negativity condition is transferred to the symmetric linear mapping 

S+A*SA--A*A--A*A = (A* S*—A* A--£~* A)+5--A* 5 _1 A 

2 2 v 2 A 2 ; 4 

for every A. Therefore, the linear mapping S — \ AS~ 1 A* has to be positive. It is 
easy to verify that S - \A* S' 1 A = |5 - \W S^W* . 

6.4. n-monotonicity of a non-symmetric 2x2 matrix. 



□ 



Proposition 6.7. Lef A be a 2 x 2 matrix with a positive definite symmetric part 
S and a nonzero skew-symmetric part 



W = 



-r 
r 



-1 

1 



= rJ. 



Let us introduce an angle 9 between and \ such that 



\r\ = VdetS tan (9 
TTzen, ^4 is n-monotone if and only if 

n8 < 7T. 

Proof. Let xi, £2, • • ■ , x n be n points in X. To study the n-monotonicity of A, 
we have to verify that the sum 



y^+i - xi,Axi) 

i=l 

is negative under the closure hypothesis x n +i = x±. Because of the linearity of A, 
this sum is invariant by translation. Let us regard x n as the origin in X, then the 
sum reduces to 

n— 1 n— 1 n— 1 



The symmetric part S of A being positive definite, we can substitute each Xi by 
S~*x;, and the sum becomes 



n-l 



n-l 



%/detS 



where we have used the identity S 2 J5 
of A is governed by the variable t 



/detS 



J. Therefore, the n-monotonicity 



/detS 



It is equivalent to study the positivity 
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of the real symmetric (2n — 2) x (2n — 2) matrix 



21 
-I-tJ 






-I + tJ 

21 






21 -I + tJ 
-I-tJ 21 



or of the Hermitian (n — 1) X (n — 1) complex matrix 



-(l + it) 






-(1-tt) 

2 









2 

-(l + it) 






-(1-it) 
2 



The eigenvalues of this tridiagonal Hermitian (ra — l)x(n— 1) matrix are well- 
known ( 11271 ). They take the (n — 1) real values 





2 






cos 




cos 6 


— cos 




cos 6* 







is ( 



indexed from A; = 1 to k = n — 1. The n-monotonicity is ensured if the smallest 
eigenvalue —^-3 [cos 6* — cos f — )] is positive, and the conclusion follows. □ 

6.5. Maximality. 

Definition 6.8. A n-monotone multifunction T is maximal if there is no proper 
n-monotone enlargement of T. 

6.6. Fitzpatrick's sequence. 

Definition 6.9. For n > 2 and (x, y) E X x Y, let (xi, yi) be n — 1 elements of 
G(T), indexed from i = 1 to n — 1, define ?/„) = (x, y), and close the loop 
by (x„+i, = (a^iij/i)- Fitzpatrick's sequence ([4 1) is defined by 

n 

FT.n{x,y) = (x,y) + sup "^(xx+i - x\,y\). 

Remark 6.10. For n — 2, we recover the function Ft.2, originally proposed by 
Fitzpatrick to study monotone operators. 

Proposition 6.11. Fitzpatrick's sequence is increasing: for n > 3, 

F T , n (x,y) > F T>n -i(x,y). 

Proof. While evaluating the supremum defining Fr n (x,y), choose a sequence 
such that x n -\ — x n -2 and compare with Ft, u -i{x, y). □ 
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Proposition 6.12. Every function of Fitzpatrick' s sequence is bounded from below 
by the duality product: 

F T , n {x,y) > i x >y) 

and equality is attained if and only ifyG Tx. 

Proof. This is a generalization of Proposition l5.10l bv Proposition l6.11l □ 

Theorem 6.13. Let T : X — > 2 Y be maximal n— monotone, for some n > 2; 
then, Ft.u represents T: 

(i) G(T) = {(x,y) eXxY\ F T , n (x,y) = (x,y)} 

(ii) outside G(T), FT <n {x,y) > (x,y) . 

Proof. We refer the reader to (E] |5]), where it is shown that Ft, u is well suited to 
study n-monotone operators. □ 

Remark 6.14. FT >n is globally lsc and convex. 

Proof. From the definition, as Ft, 2 was globally lsc and convex (Remark 15. 121 ). 
Fr,n is the upper hull of a family of continuous and affine real-valued functions 
and therefore is convex and lsc on the product space X x Y. □ 

We present next a recursive formula allowing in certain cases an iterative com- 
putation of Fitzpatrick's sequence of a monotone operator. 

7. Recursion formula 

Proposition 7.1. LetT : X — > 2 Y be a (n+l)-monotone multifunction (n > 2), 
then 

V{x, y)eXxY, F T , n+ i {x, y) = sup [F T , n (£, y) + {x - £, if)]. 

Proof. For n > 2 and (£, 77) £ X x Y, let (xi, yi) be n — 1 elements of G(T) 
indexed from i — 1 to n — 1. By definition, 

FT, n (Z>v) + ( x ~ Z>v) = 

n-2 

= SU P (z2i x i+i - Xi,yi) + (x 1} y) + (£, - x n -t,y n -i)) + (x - £,r]). 

Following the supremum over the n—1 elements (xi, yi) of G(T), by supremizing 
over (£, 77) € G(T), we recognize (JT)) 

n-2 

F T , n +i(x,y) = sup - Xj , j/j) + (x n -x n - 1 , y n - i)+(xi,y)+(x-x n ,y, 

where the role of (x n , y n ) is played by tf). □ 
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8. FlTZPATRICK'S SEQUENCE FOR A GSM 

Definition 8.1. A constitutive law associated to a multifunction T is cyclically 
monotone if it is n-monotone for every n > 2. 

Example 8.2. If X = Y is a Hilbert space, and A is a symmetric positive semi- 
definite linear mapping from X to Y, then the (single valued) multifunction T 
defined by Tx = {Ax} is cyclically monotone (see Proposition 16.5b . 

Example 8.3. If T is the subdifferential of a convex lsc potential <fi, then T = d(j) 
is cyclically monotone (|46|). 

Definition 8.4. A cyclically monotone multifunction T is maximal if T is cycli- 
cally monotone and no proper extension of T is cyclically monotone. 

8.1. Recovery of the bipotential for a GSM. A convex potential which has at 
least one finite value (not identically +oo) is called proper. When the constitutive 
law of a GSM is described by a convex, lsc, and proper potential <fr, the multi- 
function T = d<p is maximal monotone and cyclically monotone, hence maximal 
cyclically monotone. For n > 2, every n-monotonicity is captured ([4]) by Fitz- 
patrick's function FT, n (x,y). Moreover, Fitzpatrick's sequence admits a point- 
wise limit 

Ft,oo{x,v) = supF TiTl (x,y) 

n>2 

which is nothing else than the sum of the potential and its conjugate (B1l8l [T3l ) 
Ft,oo{x,v) = cp{x) + 4>* {y) 

and we recover the GSM-specific separated bipotentials. 

Conversely, when a multifunction T is maximal cyclically monotone, a con- 
structive theorem flU) due to R.T. Rockafellar (04] HU) and J.J. Moreau (ED) 
proves the existence of a convex, lsc, and proper potential <fi such that T = d<fi. 
The method for retrieving <fi(x) consists in fixing y in Ft }0 o {x, y). By duality, the 
method for retrieving <fr*(y) consists in fixing x in FT l0 o{x,y). Actually, the se- 
quence constructed by Fitzpatrick is a genial rewriting of the Moreau-Rockafellar 
theorem. 

8.2. Fitzpatrick's sequence for an indicator function. When the potential <j> is 
the indicator function %k of a convex K, then Fitzpatrick's sequence of the multi- 
function T = d<p is reduced to a single element (|4)) 

\fn > 2,F T ,„(x,y) = F T<00 {x,y) = ijc(x) +i* K (y). 

8.3. Fitzpatrick's sequence for a support function. By duality, the same equal- 
ity holds ([4|) when the potential <fi is the support function i* K of a convex set 
K. 

In the next section we will emphasize (and sometimes repeat) in the framework 
of continuum mechanics some notions and results concerning the linear monotone 
operators previously presented in sections |5]|6]and|8] 
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9. Linear constitutive laws 

9.1. Symmetric Linear Constitutive Laws. As in Example I5.8l let X and Y be 

two instances of the same Hilbert space. A linear symmetric law y — Sx models 
a Standard Material if and only if S is symmetric (i.e. coincides with its adjoint 
S*). The potential is then <f>(x) — \{x, Sx). When S is also positive definite, 
this potential is strictly convex and the symmetric linear constitutive law y = Sx 
models a well-known class of GSM: the linear elastic materials. 

Proposition 9.1. Fitzpatrick's sequence of the linear symmetric positive definite 
law y — Sx is 

F s , n {x,y) = (x,y) + 1 ^-\\S-^y-Six\\ 2 
= h x , y ) + (l-i)(h x ,Sx) + h y ,S- l y )) 

n n 2 2 

= (x, y) + (1 - -)\{y ~ Sx, S-\y - Sx)) Vn > 2 
n 2 

where Si and S~^ are the square roots of S and its inverse S^ 1 . The pointwise 
limit of Fitzpatrick's sequence is 

F s ,oo{x, y) = -{x, Sx) + - (y, S^y) . 

Proof. If in Definition 16.91 of Fs,n{ x i ll) we replace each pair (x\, y\) by the pair 
( S^x\, S~iy\\ , then we observe that it is enough to give the proof for S being 
the linear identity mapping /. But a proof for S — I was given in [4|. See also 
Propositions 1931 EOl and Example [9T5lin this Section|9] □ 

9.2. Linear Elasticity. The phenomenon described above is that of elasticity: x 
is the strain tensor, y is the stress tensor, S is the stiffness tensor, its inverse 

is the compliance tensor. The behavior of the elastic material is governed by the 
potential 4>{x) = \{x, Sx) or by the conjugate potential 4>*(y) = i (y, S~ 1 y). 

Most numerical solutions of the linear elasticity problems are obtained ( Il23ll26l 
l53ll54l[58l ) by applying variational methods: 

- the minimum principle of the potential energy, involving an integral func- 
tional of the kinematically admissible displacements fields with the inte- 
grand 4>(x) regarded as strain energy density, 

- the minimum dual principle of the complementary energy, involving an 
integral functional of the statically admissible stress fields with the inte- 
grand 4>* (y) regarded as stress energy density, 

- a primal-dual two-field variational principle, involving an integral func- 
tional with the integrand cf>(x) + 4>*(y), 

- three-field variational principles, etc. 

The bipotential approach, very close to the primal-dual two-field methods, sug- 
gests intermediate variational principles involving integral bifunctionals dl22l[39l 
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[50] [60 1) with an integrand 

F 3>n ( x >v) = ( x >y) + (!- ^M x ) + ^*(y) - ( x ,v)l 

9.3. Strictly Monotone Non-Symmetric Linear Constitutive Laws. A linear 
law y — Ax can be monotone without being symmetric, it is only necessary that 
the symmetric part S — \ (A + A*) of A is positive. It will be strictly monotone 
if S is positive definite. 

Proposition 9.2. Fitzpatrick's function of a strictly monotone linear constitutive 
law y — Ax is 

F A , 2 (x,y) = (x,y) + \{y- Ax, S'^y - Ax)). 
Proof. If yi — Axi, the infimum of 

(x - xt, y - y%) = (x - x%, y - Ax) +(x-xt,S(x- xi)) 

is attained for x\ solution of 2S(x — X\) = Ax — y. Therefore, the first element 
of Fitzpatrick's sequence is 

F A , 2 (x,y) = (x,y) + \{y - Ax, S'^y - Ax)). 

This function is a bipotential representing the non-associated linear law y = Ax. 

□ 

Proposition 9.3. Fitzpatrick's function of a (not strictly) monotone linear consti- 
tutive law y = Ax is 

Fa,2 (x, y) = (x, y) + ^ (£, (y - Ax)) + i R{S) (y - Ax) 
with £ any solution of the linear equation 2St; = y — Ax. 

Proof. If the symmetric part S of A is positive, but not positive definite, the linear 
equation 

2S(xi — x) — y — Ax 

can have solutions x\ — x = £ if and only if y — Ax belongs to the range R(S) of 
S. Therefore, the first element of Fitzpatrick's sequence is 

Fa,2{x,v) = (x,y) + -(£, (y - Ax)) +i R ^(y- Ax) 

and this value does not depend on the choice of the solution £ of the linear equation 

2S£,=y-Ax. □ 
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9.4. Strictly n-Monotone Linear Constitutive Laws. 

Proposition 9.4. Let n > 2, and Abe a strictly n-monotone linear mapping, then 

FA,n{x,y) = (x,y) + i(y - Ax,H^{y - Ax)) 

where the linear symmetric mappings Hk indexed from k ~ 2 to n are generated 
by the recursive formula Hk+i = S — \A* H^ 1 A, with an initial value H2 = S. 

Proof. The function $(.t, y) = Fa : u{x, y) — (x, y) is the supremum of 

n-2 

y^(^i+i - Xi,Axi) + (xi - x,y) + (x - x n -\,Ax n -x). 

i=l 

Let us regard x as the origin in X and set Xi — x + Zi, this translation leads to 
supremize 

n-2 

-Zi,Ax + Azi) + (zi,y) - (z n -i,Ax + Az n -i). 

i=i 

n-2 

The sum — z i, Ax) reduces to (z n -i — z\,Ax), therefore 

i=l 

n-2 

®(x,y) = sup[^(z i+ i - Zi,Azi) + (zi,y- Ax) - (z„_i, Az n -i)] 
z > i=i 

n — 2 n— 1 

= sup[}^(z i+ i,Azi) - 22(zi,Szi) + (zx,y - Ax)]. 

Zi i=l i=l 

Due to the strict monotonicity of A, the bracketed term is a strictly concave and 
differentiable function in (zi, Z2, ■ . ■ , z n -i). Then, by differentiation, the maxi- 
mum is attained for (z\, Z2, . . . , z n -i) solution of the linear system 

—25*21 + A* zi = Ax — y 

Azi-i - 2Szi + A*z i+ i = for i = 2 to n — 2 

Az n ^2 — 2S'z n _i = 0. 

When these stationarity conditions are satisfied, the quadratic part of the above 
concave function is the additive inverse of the half of the linear part, and ll57l 

Fa,u (x, y) = (x, y) + - (21, y - Ax) . 

It remains to express z\ in terms of y ~ Ax. Let us rewrite the above linear system 
as 

Szi - \a*z 2 = i (y - Ax) 
Szi — ^A* z i+ i — -Az.i-i for i = 2 to n — 2 

Sz n ~i = \Az n ^2- 
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Clearly, this linear system can be solved by backward substitution. Let us ini- 
tialize by H2 = S a sequence Hk (indexed from k = 2 to n) of invertible linear 
symmetric mappings, and set 

then successively: 

H„_ l+1 = S - -A*H~ 1 -A for i = n - 2 to 1 

4 

We end with 

H n zi = ^(y- Ax) 

and the conclusion follows. □ 

Example 9.5. If A is the identity /, the sequence Hk is generated by the recursive 
formula Hk+i — I — \H^ . Because of the initial value H2 — /, each Hk is 
spheric, and we can set Hk = ajtJ. The sequence of scalars au is generated by 
the recursive formula ak+i = 1 — 3^-. The homographic function 1 — admits 
the unique fixed point a — 5. As a property of the complex projective line P 1 (C), 
the sequence fik = — is an arithmetic progression. The common difference is 
equal to 2: fik+i = ^fc + 2. The initial term being 012 = 1, i.e. $2 = 2, the solution 
is Pk = 2(fc — 1) i.e. otk = \ -^zj ■ Thus H^ 1 = 22^7 and we retrieve the result 
given in [4] 

Fi, n (x,y) = (x,y) + (1 - -)~{y- x,y- x). 

n 2 

Example 9.6. If A = S, the sequence Hk is generated by the recursive formula 
Hk+i — S — jSH^ 1 S. Because of the initialization H 2 = S, each Hk is pro- 
portional to S, and we can set Hk = akS. The sequence of scalars ctk is as in 
Example 19.51 Thus H~ l = and we retrieve the result of Proposition 

eh 

Fs, n (x,y) = (x,y) + (1 )-{y - Sx,S^ 1 (y - Sx)). 

n 2 

Remark 9.7. The symmetric linear mappings Hk are positive definite. 

Remark 9.8. The recursive construction Hk+i — S — \A* H^ 1 A stops if Hk 
ceases to be positive definite, i.e. if A fails to be (fc + l)-monotone. 

Remark 9.9. The function Fa, 3 captures the 3-monotonicity (Example 16.4b of 
the linear mapping A. 

Remark 9.10. The special recursion formula H n+ \ = S — \ A* H~ x A must be in 
accordance with the general recursion formula ( Proposition |7J} 

FA,n+i{x,y) = sup [F A ,n{£,,y) + (x-£,A£,)] . 
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To evaluate this supremum, let us regard x as the origin in X and set £ = x + (. 
This translation leads to supremize the function 

(x + (,y)- (C, Ax + A() + 1 (y - Ax - AC, H~\y -Ax- AC)) 
with respect to £. The quadratic part of this function reads 

- (CSO + ^ACH-'AC) = -{C,H n+1 C). 
The linear part of this function reads (£, 0) with 

z = 1J Ax \A*H-\y Ax) = (I - ^A'H^Xy Ax) 
The constant part of this function reads 

(x, y) + \{y- Ax, H-\y - Ax)). 

If A is a strictly (n + l)-monotone linear mapping, the function to maximize is 
strictly concave. The maximum is attained for C solution of the linear system 

2-ffn+lC = z - 

When this stationarity condition is satisfied, the quadratic part is the additive in- 
verse of the half of the linear part, and 

F Ain+1 (x,y) = (x,y) + j(y-Ax,H- 1 (y-Ax)) + j(z,H-l 1 z) 
has to coincide with 

F A . n+1 {x, y) = (x, y) + ^{y- Ax, H'l^y - Ax)). 
Therefore, the sequence H^ 1 follows the formula 

H- 1 + ^H-'A^-^il \a*H- x ) = H-l, 

which allows us to compute recursively K n+ \ = H~+ x from K n — H~ x (hence 
from K 2 = S^ 1 ) by solving the equation 

K = K n + (I- X -K n A)K(l - \A*K n ) 
using a fixed point algorithm initialized by K n . 

10. Study of linear monotone coaxial constitutive laws 
10.1. Coaxial constitutive laws. 

Definition 10.1. Let X and Y be the 6-dimensional Euclidean space of real sym- 
metric 3x3 matrices (with e as identity matrix and (x,y) = tr(xy) as duality 
product). The variables x and y can be regarded as strain and stress tensors. A 
constitutive law relating x and y is coaxial if x and y have the same eigenvectors. 
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Proposition 10.2. Under the additional assumption of linearity, the coaxial con- 
stitutive laws have the general form 

y = [tr(kx)\ e + 2 fix 

where fi is a scalar and k is a symmetric 3x3 matrix. 

Proof. Let us begin by making two remarks 

(i) if x = e, then y is spheric, 

(ii) for any unit vector u, if x = uu*, then y = a(u)uu* + (i(u)e where a 
and j3 are two scalar functions. 

Consider 3 unit vectors u, v, w constituting an orthonormal basis. Taking 

x = uu* + vv* + ww* = e, 
we can conclude from the linearity assumption that 

a{u)uu* + a(v)vv* + a(w)ww* + [(3(u) + j3(v) + j3(w))e 

is spheric. This implies that the scalar function a has a constant value that we will 
denoted by 2/i. Applying now the linearity assumption to the generic case 

x = \\uu* + X 2 vv* + X 3 ww*, 

we conclude that y = 2[ix + -f(x)e where 7 is a scalar linear function. But any 
linear scalar function of x can be expressed as tr(kx) where A; is a symmetric 3x3 
matrix, and the conclusion follows. □ 

Remark 10.3. A linear coaxial constitutive law involves 7 coefficients: the scalar 
fi and 6 independent coefficients of the symmetric matrix k. 

Remark 10.4. A linear coaxial constitutive law is not symmetric except when the 
matrix k is spheric (k = Ae), in which case A and fi are Lame's coefficients of 
Hooke's elastic constitutive law 

y = \(trx)e + 2/j.x. 

It is a well-known result that Hooke's constitutive law is positive if and only if 

3A + > and fj, > 0. 

10.2. Symmetric part of a linear coaxial constitutive law. If the deviatoric part 
h of the matrix k is not 0, set 

fc = Xe + h with A = - (tr fc) and tr h = 0. 
o 

The symmetric part of the linear coaxial constitutive law is then 

Sx — i [tr(fcx)] e + i(tr.x)fc + 2[ix 

= A(trx)e + 2[ix + ^-(trx)h + ^ [tr(/ix)] e. 
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10.3. Monotonicity of linear coaxial constitutive laws. 

Proposition 10.5. A linear coaxial constitutive laws is monotone if and only if 

(131132) 

g 

3A + 2/i > 0, fi > and tr (h 2 ) < -/i(3A + 2/i). 

Remark 10.6. In addition to the usual conditions on Lame's coefficients, the 
mono- tonicity condition of the coaxial constitutive law demands that the devi- 
atoric part h of the matrix k has not to be too large. 



Proof. If the deviatoric tensor h is not vanishing, let us express the matrix of the 
symmetric part S of the linear coaxial constitutive law in an orthonormal basis 
constituted of four deviatoric tensors d\, d%, d$, d± orthogonal to h, d<^ = and 



d R = 



V3' 



For i = 1 to 4, Sdi = 2jj,di, otherwise Sd§ = 2/id5 



and 



Se — (3 A + 2/i)e + The matrix of S is therefore the symmetric block matrix 



2(ih 













2 

3A + 2/i 



2 M 



where I4 is the 4x4 identity matrix. The 6x6 symmetrical matrix representing 
S is positive if and only if the scalar /t and the 2x2 symmetrical matrix 



2/i 



3A 



2// 



are positive. The additional conditions for s to be positive are 



3A + 2/i > and -^\\h\\ 2 < 2/i(3A + 2/i). 



□ 



10.4. A characteristic angle of monotone linear coaxial laws. Let us rewrite 
the inequality limiting the magnitude of h as follows: 

\\h\\ < -^V2/i(3A + 2/i). 

We are led to introduce an angle 9 between and ^ such that 

\\h\\ = -^= n /2m(3A + 2/x) sin 6. 
v3 

The smaller 9 is, the closer is the linear coaxial law to Hooke's constitutive law. If 
29 < 7r, then the law is strictly monotone. 
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10.5. Condition of n-monotonicity for a linear coaxial law. 

Proposition 10.7. A linear coaxial law is n-monotone if and only if 

3A + 2u > 0, fx > and n9 < tt. 

Proof. A monotone linear coaxial law is represented by a 6 x 6 block matrix, the 
diagonal of which is compound of a positive definite spheric 4x4 matrix 2/x/4 
and a 2 x 2 matrix 



2[i 





3A + 2/i 



-1 



s + rJ 



with 



2 



— \\h\\ = Vdet s tan 6. 



According to Proposition 16.71 the additional condition for the constitutive law to 
be n-monotone is n9 < tt. □ 

10.6. Condition of cyclic-monotonicity for a linear coaxial law. 

Proposition 10.8. A linear coaxial law which is cyclically monotone reduces to 
Hooke 's constitutive law. 

Proof. If a linear coaxial law is cyclically monotone, then n9 < tt for every integer 
n larger than 2. The field of the real numbers being Archimedean, this is possible 
only for 9 = 0, i.e. h = 0. Therefore the linear coaxial law reduces to Hooke's 
constitutive law. □ 

10.7. Fitzpatrick's sequence of a strictly n-monotone linear coaxial law. 

Proposition 10.9. Let A be a strictly n-monotone linear coaxial law, then Fitz- 
patrick's sequence of A is finite and constituted of the n — 1 functions Fa,h defined 
for k — 2 to n by 



F A>k {x, y) = tr(xy) + -tr [(y - Ax)H k 1 {y - Ax)} 



with 



2/iQ!fc/4 






h k _ 




n i sin(ke) i „ 

u 2 sin((fc-l)0) cos0 _ 



Proof. The orthonormal basis d\, di, d%, d±, = j^, rf 6 = already chosen 
for representing S (Proposition ll0.5l i. leads to represent the linear mapping A by 
the 6 x 6 matrix (Proposition 1 10.71 ) 



A = 



2al 4 




where a is the 2x2 matrix 



2u 

y/3\\h\\ 





3X + 2fi 
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To express the k th Fitzpatrick's function we have to solve the recurrence equa- 
tion Hk+i = S — jA* H^ 1 A. We can transform this problem into two smaller 
problems: 

(i) find the solution for the 4x4 spherical matrix 2^1^, 



(ii) find the solution for the 2x2 matrix a. 

Up to the multiplicative factor 2[i, the first problem was already solved in Example 
19.51 the solution is 2iia.kli = A* fc~T^ 4- ^ remains to solve the recurrence equation 
hk+i — s — ja*h~ k 1 a initialised with h% = s. By induction, we observe that hk is 
proportional to s. Let us set hk = \lkS- Therefore, 7fc is satisfying the recurrence 
equation 7fc+i = 2 — co ^ g with 72 = 2. As one can easily check by induction, 

the solution is 71. = ■ B ^ k9 X a \ —^s, which achieves the proof. □ 

sin((/c— 1)0) cos 9 ' * 

Remark 10.10. If we introduce the variable X — cos 0, then X~/k = ul llx] ' s 
a quotient of consecutive Chebyshev polynomials of the second kind. The above 
recurrence equation satisfied by the 7^ and the initial value 72 = 2 follow from 
the recurrence relation and initial conditions 

U {X) = 1, Ui(X) - 2X, U k (X) = 2XU k - 1 {X) - U k - 2 {X) 

defining these polynomials. The well-known solution of this recurrence problem 
is Uk{X) = "^f " 1 "^^ , and we recover the above trigonometric expression of 7/.. 

Remark 10.11. As approaches zero, jk approaches = 2a k , and hk ap- 
~2fj, 



proaches a k 



3A + 2/i 



in accordance with Example 



11. Revising elasticity theory-Return to Hooke 

11.1. What did Robert Hooke say? First, he noticed that the extension of a 
spring is proportional to the weight hanging on it. In 1660, from this experimental 
observation, he modeled the behavior of elastic materials by stating the law "UT 
TENSIO SIC VIS" (Latin: as the extension, so the force). He published it in the 
anagram form "ceiiinosssttuv", whose solution he gave in 1678. 

Nowadays, this law is interpreted as stating that the stress tensor y is an isotropic 
linear function of the strain tensor x: 

y = X(tix)e + 2 fix 

where A and /1 are coefficients (Lame's coefficients). 



In the next subsection, we discuss this interpretation of Hooke's prescription. 
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11.2. Interpretation of Hooke's prescription. Hooke's prescription was one- 
dimensional. At that time, to model the behavior of elastic bodies, Robert Hooke 
did not have at his disposal the mathematical concepts of vectors and tensors. 
How to interpret his prescription as a 6-dimensional one between strain and stress 
tensors? From Hooke's "UT / SIC ", we will retain the hypotheses of linearity 
and monotonicity. Why to add the condition of isotropy? By analogy with the 
behavior of springs, Robert Hooke claims that a force applied to an elastic body 
in some direction generates a proportional deformation in the same direction. In 
the framework of tensorial calculus, we will understand it as: if a direction is 
an eigenvector of the stress tensor, then it is an eigenvector of the strain tensor. 
Finally, we will interpret Hooke's prescription by requiring the conditions of 

- linearity, 

- monotonicity, 

- and coaxiality, 

for the relation between the strain tensor x and the stress tensor y. According to 
Proposition 110.21 these three hypotheses lead to the constitutive law: 



where A and fi are scalar coefficients (Lame's coefficients) and h is a deviatoric 
tensor. According to Proposition [TOTS] the monotonicity is insured by the inequal- 



1 1.3. Four reasons for which the coaxial law reduces to the isotropic law. 

11.3.1. Classical argument: isotropy. The isotropy assumption asks, for every 
rotation R, that the change of x in R~ 1 xR is followed by the change of y in 
R _1 yR. Therefore, the deviatoric tensor h has to satisfy R~ 1 hR = h, and is 
forced to vanish. 

11.3.2. Hermann von Helmholtz's argument: existence of a strain energy. In the 
linear case, a strain energy density exists if an only if the constitutive law is sym- 
metrical. Therefore, the deviatoric tensor h is forced to vanish. 

11.3.3. Lars Onsager's argument: symmetry. As in Remark ll0.4l the symmetry 
of the constitutive law implies that the deviatoric tensor h vanishes. 

11.3.4. Jean-Jacques Moreau's argument: cyclic monotonicity. The cyclic mono- 
tonicity leads to the existence of a strain energy, and therefore to the vanishing of 
the deviatoric tensor h (see also Proposition ll0.8b . 

Certainly, the reader is now aware that in our opinion, the most relevant argu- 
ment is the fourth (cyclic monotonicity) and not the first one (isotropy). 



y = A(tra;)e + 2fix + tr(hx)e 



ities: 
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1 1 .4. Revisiting isotropic elastic materials. In light of the previous results, we 
do believe that elasticity theory must be revisited. The constitutive laws of the 
so-called isotropic linear elastic materials have to be extended to monotone linear 
coaxial laws. The cyclic monotonicity is a too strong assumption, fc-monotonicity 
up to a finite integer n has to be considered. We do believe that this maximal inte- 
ger n is a relevant parameter characterizing the material, and we will say that the 
material is n-monotone. Classical elasticity corresponds to very large integers n 
(therefore to very small angles 9, cf !10.5b - The evaluation of the seven coefficients 
of the model must be performed by modern identification methods, allowing an 
optimal exploitation of the electronic recording of the measures ((2] 0). Forgo- 
ing the cyclic monotonicity, the existence of an elastic potential <f>{x) is lost. But 
numerical methods based on primal-dual two-field variational principles can be 
easily generalized. Simply, the integrand <fi(x) + <f>*(y) of the classical integral 
functional has to be replaced (G2[39]|50]|60)) by the bipotential 

F A>n {x, y) = tr(z, y) + -tr[(y - Ax)^ 1 {y - Ax)} 



with 







in accordance with the notation 





1 sin(n9) 1 

2 sin((n-l)0) cos 9 u ' 



2fi 



VI 

2 

3X + 2fi 



and with the definition of the angle 9 



v/3, 



(3A + 2/z) sin 9 = -y ||/i|| 



12. Conclusion 

The class of n-monotone materials for which the constitutive law is described 
by a n-monotone operator is larger than the class of Generalized Standard Mate- 
rials (GSM). The integer n can be regarded as a characteristic of these materials. 
Fitzpatrick's functions allow to construct algorithms based on primal-dual two- 
field variational principles, as easy as the GSM-specific algorithms ( Il22l[39l|50l 
60 1). To apply this modeling in fluid mechanics, it is necessary to replace the 
tensor x by the strain rate tensor and the tensor y by the Cauchy stress tensor, aug- 
mented by the spheric tensor pe, where p is the pressure. To apply this modeling 
in thermal engineering, it is necessary to replace the tensor x by the opposite tem- 
perature gradient and the tensor y by the heat flux vector. To apply this modeling 
in electromagnetism, it is necessary ([30, 3T|) to replace the tensor x by 

E 1 E 2 E 3 ' 

E 1 —B 3 B 2 

E 2 B 3 -B 1 

E 3 —B 2 B 1 
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and the tensor y by 

"0 D 1 D 2 D 3 

D 1 -H 3 H 2 

D 2 H 3 -H 1 

D 3 -H 2 H 1 

where E, B, D, H, stand respectively for the electric field strength, the magnetic 
induction field, the electric displacement and the magnetic field intensity. 

However, in mechanical and civil engineering, the constitutive laws of many 
materials (ductile metals, metal matrix composites, wet clays, plastic soils, gran- 
ular materials, etc) are not monotone. The class of Implicit Standard Material 
(ISM) is larger than the class of n-monotone materials. Every Fitzpatrick's func- 
tion is globally lsc and convex, but the modeling of ISM only requires bipotentials 
which are partially lsc and convex. It would be interesting to generalize the con- 
cept of Fitzpatrick's sequences to the case of non-monotone operators. It is also 
expected that the recurrence formula of Proposition l7.1 I passes this extension with 
few changes. 

As a pioneering attempt for modeling Coulomb's dry friction, let us consider 
the following constitutive law in an Euclidean linear space: two vectors x and y 
have the same orientation. This constitutive law is not monotone, Fitzpatrick's 
method cannot be directly applied. Does this constitutive law model an IMS? Can 
Fitzpatrick's sequence be generalized for representing it by bipotentials? If we 
only ask for a local supremum in the Definition 15.91 of Fitzpatrick's function, we 
obtain (|56|) the extremal value 

h(x,y) = ^(x,y) + ^|M||M|. 

This function is partially lsc and convex. Thanks to the Cauchy-Schwarz-Bunia- 
kovsky inequality, it is a bipotential. If ip is an angle chosen between and ir such 
that (x, y) = \\x\\ \\y\\ cos ip, we observe that 

( tp 

b 2 {x,y) = \\x\\\\y\\ I cos - 

A similar weakening in the definitions of Fitzpatrick's functions FT in {x, y) leads 
( ll56ll ) to the following increasing sequence of bipotentials 

( i> 

b n (x,y) = \\x\\\\y\\ I cos - 

The pointwise limit 

b(x,y) = \\x\\\\y\\ 

is a wellknown bipotential (|55|) representing the constitutive law asserting that 
two vectors x and y admit the same orientation; we will call it "Cauchy-Schwarz- 
Buniakovsky bipotential". 

We hope that this kind of extension of Fitzpatrick's sequences will prove to be 
very helpful to produce relevant bipotentials for representing the non-associated 
constitutive laws of the ISM evoked in subsection l4.3l 
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Finally, it is worth mentioning that the notion of bipotential can provide ad- 
ditional modeling to rate-independent phenomena with hysteretic behavior : lin- 
earized plasticity with hardening, finite-strain elastoplasticity with nonlinear and 
non-associated kinematic hardening rules, damage in ductile materials, phase trans- 
formations in shape memory alloys, delamination, ferromagnetism, superconduc- 
tivity ( fl4l[l5ll38l f43'l). quasistatic evolution of fractures, and crack propagation 
in brittle materials. 
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